EQUIVALENCE OF WEAK FORMULATIONS OF THE STEADY WATER 

WAVES EQUATIONS 



EUGEN VARVARUCA AND ARGHIR ZARNESCU 



Abstract. We prove the equivalence of three weak formulations of the steady water waves 
equations, namely the velocity formulation, the stream function formulation, and the Dubreil- 
Jacotin formulation, under weak Holder regularity assumptions on their solutions. 



1. Introduction 

We consider the classical problem of travelling waves at the free surface of a two-dimensional 
inviscid, incompressible, heavy fluid over a flat bed. Over the past few years this problem 
has attracted considerable interest, see [10] for a survey of recent developments. In the most 
direct mathematical description, the problem is to find steady solutions of the incompressible 
Euler system for the velocity field (m, v) and the pressure field P, together with relevant 
boundary conditions, in an unknown domain in the plane. In this form the problem is difficult 
to treat mathematically, and other, more convenient reformulations, have been used in the 
literature. One such reformulation involves a stream function ip, whose existence is ensured 
by the incompressibility condition. Namely, the stream function ip satisfies a semilinear 
elliptic equation (since vorticity may be present in the flow), together with suitable boundary 
conditions. This reformulation still retains the difficulty of being a free-boundary problem. 
The most general approach to overcoming this difficulty uses a change of variables introduced 
by Dubreil-Jacotin in 1934, which transforms the problem to an equivalent one in a strip. 
The Dubreil-Jacotin formulation is a cornerstone of the large and growing literature on large- 
amplitude travelling water waves with vorticity that began with the work of Constantin and 
Strauss [2]. (For other methods, of more restricted applicability, of transforming the free- 
boundary problem for the stream function into a problem over a fixed domain, see [1], [6], 

As we discuss in Section 2, following essentially [2], the three formulations of the steady 
water waves equations just described, to which we refer to as the velocity formulation, the 
stream function formulation and the Dubreil-Jacotin (or height) formulation, are equivalent 
when considered in the classical sense. However, the recent study [3] by Constantin and 
Strauss has raised the question whether the three formulations are still equivalent when con- 
sidered in a suitable weak sense which we describe in Section 3. The main result of [3] is a 
global bifurcation theory for solutions, periodic in the horizontal direction, of class C^'" for 
some a G (0, 1) of the weak Dubreil-Jacotin equation. Such solutions would formally corre- 
spond to solutions {u,v) and P, also periodic in the horizontal direction, and of class C°'°, of 
the weak velocity formulation, in a domain with C^'" boundary. However, the key question 
of the equivalence of the two weak formulations under these regularity assumptions is left as 
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an open problem by [3]. (A result on the equivalence of the formulations is given in [3], but 
under different, and not the most natural, regularity assumptions.) 

The main result of the present paper. Theorem 1, which is given in Section 4, provides 
an affirmative answer to the above open problem, albeit only in the case when the Holder 
exponent satisfies a > 1/3. More precisely. Theorem 1 proves that, under the above regularity 
assumptions, the weak velocity and the weak stream function formulations are equivalent for 
any a G (1/3,1], while the weak stream function and the weak height formulations are 
equivalent for any a G (0, 1]. An important consequence of our result is that, at least in the 
case a > 1/3, the solutions constructed in [3] of the Dubreil-Jacotin formulation are relevant 
(in the sense that they give rise to corresponding solutions) for the velocity formulation of 
the steady water waves equations. Our result is in the same spirit as, and its proof is inspired 
by, the Onsager conjecture as proved (partially) in [5]. The Onsager conjecture is, essentially, 
the statement that solutions of the time-dependent incompressible Euler equations on a fixed 
domain (in dimension three, with no external forces), of class C"'" in the space variables for 
each value of the time variable, conserve their energy in time if a > 1/3 and may fail to do 
so if a < 1/3. The paper [5] proves that a > 1/3 implies conservation of energy (and leaves 
open the reverse statement in the conjecture). As in [5], our proof is based on regularizing 
the equations and, roughly speaking, the assumption a > 1/3 is used in an essential way to 
show that certain remainder terms converge to as the regularization parameter tends to 0. 
An important problem left open by the present paper is that of whether the weak velocity 
formulation and the weak stream function formulation are also equivalent in the case when 
the Holder exponent satisfies a < 1/3. 

2. Classical formulations of the steady water waves problem 

2.1. The velocity formulation. We consider a wave travelling with constant speed and 
without change of shape on the free surface of a two-dimensional inviscid, incompressible 
fluid of unit density, acted on by gravity, over a flat, horizontal, impermeable bed. This 
means that, in a frame of reference moving at the speed c of the wave, the fluid is in steady 
flow in a fixed domain. Let the free surface be given hj y = ri{x), for some function 77 : M — ?■ M, 
and the fiat bottom be given by ?/ = 0, so that the fluid domain is 

= {{x,y) : X eR,0 <y < r]{x)}. 

Then the fluid motion is described, see ^ for details, by the following equations and boundary 
conditions for a steady velocity field {u,v) and a pressure field P in D^: 



(2.1a) Ux + Vy = Drf, 

(2.1b) {u — c)ux + vUy = —Px in D^, 

(2.1c) {u — c)vx + vvy = —Py — g in D^,, 

(2.1d) ^; = ony = 0, 

(2.1e) V = {u — c)rix ony = ri{x), 

(2.1f) P = Patm ony = 7]{x), 



where Patm is the constant atmospheric pressure and g is the gravitational constant of acceler- 
ation. More precisely, fl2.1al) describes mass conservation, (]2.1bl) - fl2.1cp describe momentum 
conservation, (12.1dp expresses the fact that the velocity at the bottom is horizontal, (I2.1ep is 
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the kinematic condition that the same particles always form the free surface, while fl2.1f|) is 



the dynamic condition that at the free surface the pressure in the fluid equals the constant 
atmospheric pressure. This is a free-boundary problem, because the domain is not known 
a priori. The system (12.11) will be referred to as the velocity formulation of the steady water 
waves equations. Throughout the paper we assume that, in the moving frame, the horizontal 
velocity of all the particles is in the same direction. For definiteness, we assume that 

(2.2) u < c in Dr,. 

(All the results discussed in the paper have corresponding analogues if instead of (12. 2 p one 
assumes that u > c in D^.) 

For the remainder of this section we describe informally, following [2] , two other equivalent 
formulations of (12. ip . assuming that the solutions are smooth enough. The equivalence of 
these formulations under weak regularity assumptions is the main aim of the paper, which 
will be addressed in the subsequent sections. 

2.2. The stream function formulation. Suppose that (12. ip and (12. 2p hold. Equation 
(12. lap imphes the existence of a function ijj in D^, called a (relative) stream function, such 
that 

(2.3) il^y = u — c, tpx = —V in D^. 

The boundary conditions ( ]2.1dp and (I2.1ep imply that ip is a. constant on each of ?/ = and 
y = r]{x). Since ip is only determined up to an additive constant, one can assume that ip = 
on y = ri{x), and then we obtain that there exists a constant po such that ip = —po on y = 0. 
The condition (12. 2 p can be rewritten as 

(2.4) ipy<0 in D^, 

a consequence of which is that po < 0. After expressing the left-hand side in fl2.1bp and (I2.1cp 



in terms of ip, differentiation of the first of these equations with respect to y and of the second 
with respect to x allows us to eliminate the pressure, leading to 

(2.5) {AiP),iPy = {AiP)yiP, in D^, 

where A denotes the Laplace operator. Note that (12. 4p shows that all the level sets of ip are 
graphs over the x coordinate, and (12. 5p then implies that Aip is constant on each level set of 
Ip. Thus there exists a function 7 : [0, —po] — t- M such that 

—Alp = 'j{ip) in D^. 

(Since the quantity 00 '= — Uy = —Aip has a physical interpretation as the vorticity of the 
flow, the function 7 is customarily referred to in the literature as the vorticity function.) Let 

(2.6) r{p) = / 7(-s) ds for all p G [po, 0]. 

Jo 



It is then easy to verify, using fl2.1bp - (l2.1cl) . that 



(2.7) P +^\Vip\'^ + gy -T{-ip) = constant in D^. 

In view of fl2lf|l and the fact that = on y = ri{x), it follows that 

\Vip\'^ + 2gy = Q ony = 7]{x), 
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for some constant Q. We have therefore obtained the stream function formulation of the 
steady water waves equations, which is to find a domain and a function %Ij in Dj^ such that 



(2.8a) A?/- = -7(V') in D^, 

(2.8b) ip = —po on y = 0, 

(2.8c) ip = on y = ri{x), 

(2.8d) \V^\^ + 2gy = Q on y = T]{x), 

for some constants Pq < and Q, and some function 7 : [0, —po] — > 

Conversely, suppose that ip satisfies (12. 8 p and (12. 4p in a domain D^. Then one can define 
in a vefocity field {u, v) by (12. 3p and a pressure field P by (12. 7p with a suitable choice of 
the constant in the right-hand side, and easily check that (12. ip and (12. 2p hold. 

2.3. The height (or Dubreil-Jacotin) formulation. An elegant way to overcome the 
difficulty that in (12. 8 p the fiuid domain needs to be found as part of the solution has been 
first observed by Dubreil-Jacotin: the fact that is constant on the top and the bottom of 
Drj can be used to transform (12. 8 p into a nonlinear elliptic boundary-value problem in a fixed 
domain. More precisely, suppose that (12. 8 p and (12. 4 p hold, and let us consider the partial 
hodograph (or semi-Lagrangian) mapping 

(2.9) (x, y) (g,p) = (x, -^p{x, y)), 

which is, as a consequence of (12. 4p . a bijection between and the closure of the strip 

R = {iq,p) ■.qeR,peipo,0)}. 
Then the inverse mapping, from R to -D^, necessarily has the form 

(2.10) {q,p) ^ {x,y) = {q,h{q,p)), 

for some function /i : — M. More precisely, the following two relations hold: 

(2.11) — ip{q, h{q,p)) = p for all {q,p) G R, h{x, —ip{x, y)) = y for all (x, y) G Drj. 

(These relations show that, for each {q,p) G R, one may interpret h{q,p) as the height of the 
streamline ip = —p above the point (q', 0) on the bed.) The condition (12.41) can be expressed 

as 

(2.12) /ip > inR. 
Note also that 



(2-13) hq = — ^, hp = — — , = ijjy 



hp hp 



and 



1 hg^ , 1 

P 

hp--'' - hp'^'' ^"-^^^ "V^/^'^^- 



(2.14) = dq- ^dp, dy = -^dp, dq = - ^dy, dp = — ^9^. 
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Using these identities, one can easily reformulate fl2.8p as the following system for the function 
h defined above: 

(2.15a) (1 + hl)hpp - 2hghphgp + /ij/ig^ = -li-p)hl in R, 

(2.15b) h = on p = po, 

(2.15c) l + hl + {2gh-Q)hl = on p = 0. 

This is the height (or Duhreil-Jacotin) formulation of the steady water waves equations. 

Conversely, suppose that h satisfies (12.151) and (12.121) . Let r/ : M — )■ M be given by 77(g) = 
h{q,0) for all g G M. Then (I2.12p implies that (g,p) ^ {x,y) = {q,h{q,p)) is a bijection 
between R and Drj. Defining ip by (12. lip , the formulae (I2.13p -( l2.14p are valid, and one can 
easily deduce from (12J[5|) and ( 12T2|) that ([21]) and ([MD hold. 



3. Weak formulations of the steady water waves problem 

3.1. Weak velocity formulation. For sufficiently smooth functions rj, u,v, and P, (12. ip is 

easily seen to be equivalent to 

(3.1a) {u — c)x + Vy = in D^, 

(3.1b) ^^u-cy), + {{u-c)v)y = -P, inD„ 

(3.1c) {{u- c)v)x + {v'^)y = -Py - g in D^, 

(3. Id) v = ony = 0, 

(3.1e) V = {u — c)r]x ony = ri{x), 

(3.1f) P = Patm ony = 7]{X). 

However, (13. ip may be given a meaning for functions of weaker regularity than those of (12. ip . 
namely by interpreting (I3.1ap -( l3.1cp in the sense of distributions. Of particular interest for 
us wi ll be solutions of dtl]) with r] G Ci'"(M) and (m,w,P) G C°'"(A^) for some a G (0,1], 
with (13.1dp -( l3TTfj) being satisfied in the classical sense, and (I3.1ap -( 13.1cp being satisfied in the 



sense of distributions. (Under the same regularity assumptions, it is not clear how to give a 
meaning directly to (12. ip . because the multiplication of a distribution by a function of finite 
differentiability is not well defined.) 

3.2. Weak stream function formulation. For sufficiently smooth functions ip and 7, the 
algebraic identity 

ii^.i^y). ^l)y - (r(-V^)), = i^yiM + im 

shows that, in the presence of (12. 4p . (12. 8p is equivalent to 

(3.2a) (^,^,), - ]^{iPl - ,pl)y - (r(-^)), = in D„ 

(3.2b) = -pq ony = 0, 

(3.2c) Ip = on y = r]{x), 

(3.2d) \Vip\'^ + 2gy = Q on y = r/(x). 
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Again, fl3.2ap may be required to hold in the sense of distributions. We will be interested in 
soluti ons o f O with r] e Ci'"(M), ip G and T G CO.«([po,0]) for some a G (0,1], 

with f l3.2bp -( r3.2dp being satisfied in the classical sense, and fl3.2ap being satisfied in the sense 



of distributions (with ipxyi^y being understood in the classical sense). 

3.3. Weak height formulation. For sufficiently smooth functions h and 7, the algebraic 
identity 

shows that, in the presence of fl2.12p . (12.151) is equivalent to 
(3.3a) l-^-^ + m] +IM =0 ini?. 



, 2hl '"^^ J [h„j 
(3.3b) h = on p = pq, 

(3.3c) __l + ^/,_| = onp = 0. 

We will be interested in solutions of (13. 3p with h G C^'"(i?) and F G C°'"([po; 0]) for some 
a G (0, 1], with f l3.3bp - (l3.3cp being satisfied in the classical sense, and (I3.3ap being satisfied 



in the sense of distributions (with hp, hq understood in the classical sense). 



4. The equivalence of the weak formulations 

Weak solutions (in the sense described in the previous section) of the steady water waves 
problem have been studied only very recently in [3]. That paper deals with waves which 
are periodic in the horizontal direction, the subscript per being used in what follows to 
indicate this periodicity requirement. In |3] the authors develop a global bifurcation theory 
for weak solutions of (13. 3p with h G Cp'^{R), under the assumption F G C°'"([po5 0]), for some 
a G (0,1). These would formally correspond to solutions of the weak velocity formulation 
with r] G Cp'^iM.) and u,v,P G C^'^{D,^). However, no rigorous proof of this equivalence 
is given in [5]. The only result there on the equivalence of the weak formulations, see [21 
Theorem 2], is the following: 

Let < a < 1 and r = . Then the following are equivalent: 

(i) the weak velocity formulation 113. 1\) together with ^2.^) . for rj G Cp;,"(]R) and u,v,P E 
W^;,:{D,) C COi?(A;); 

(ii) the stream function formulation /i2.8\) together with \2.4\) , for 7 G L''[0, — po]? V £ 
ClCM) and ^ G Wl^D^) C C1^?(A;); 

(iii) the weak height formulation 113. 3\} together with 112. 12\} . for F G iy"^'^[po;0] and h G 

As one can see, in the above result the velocity field {u,v), the pressure P, the stream 
function ip, the height h, and the function F, are assumed to have more regularity, namely an 
additional weak (Sobolev space) derivative, than one would really like. 
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Our main result, given below, proves the equivalence of the weak formulations under the 
'right' regularity assumptions, albeit only for the case when the Holder exponent satisfies 
a G (1/3, 1]. (In particular, under our assumptions, the function F need not have a (weak) 
derivative.) While the weak stream function and the weak height formulations will be seen 
to be, in fact, equivalent for any a G (0, 1], it remains an open problem whether the weak 
velocity and the weak stream function formulations are equivalent for a G (0, 1/3] also. For 
simplicity, we state our result for solutions which are periodic in the horizontal direction, 
though this assumption is not essential, and the result can be easily extended to cover other 
important situations, for example that of solitary waves [7]. 

Theorem 1. Let a G (1/3, 1]. Then the following are equivalent: 

(i) the weak velocity formulation ^3.1\) together with 112. for r] G Cp;°(]R) and u,v,P E 

(ii) the weak stream function formulation h3. ^)) together with \2.Ji\} , for F G C°'"([po5 0]), 
r]eClC,m andijeCl':i(D~,); 

(iii) the weak height formulation 113. 3\} together with 112. 12\} . for F G C°'°([j9o, 0]) and h G 

A key ingredient in our proof is regularization of the relevant equations. Therefore, we 
start with some background results on regularization, valid in any number d of dimensions, 
where we use x, y and z to denote points in M'^. Let g G C^(M'^) be a given function, such 
that 

g > in W^, supp g C i?i(0), g{x) = g{—x) for all x G M'^, and / g{x) dx = 1 

and let us denote g^ == ^^i(|). Let V be an open set in W^, and consider, for any e > 0, the 

set V {x eV : dist(x, M'' \ V) > e}. For any / G Lj^^i^) and any e > such that V is 
not empty, consider in the function 

(4.1) fe^a:)"^ f*g%x) 



Q''{x-y)f{y) dy 

I g(z)f(x — ez) dz for all x G V^. 

Jb,(o) 



'Bi(O) 

For any f,g & Lf^^{V), we also introduce in the function 

(4.2) r%f,g){x)= g{z){f{x - ez) - f{x)){g{x - ez) - g{x)) dz for all x G 

^Bi(O) 

We further denote 

(4.3) R^^f^g)''^r^i^f,g)-{f-f^){g-g^) in 
Then one can easily check that we have, at every point in V^, 

(4.4) U9Y = r9' + R'U.9)- 

(It may be worth pointing out that, while it is not immediately clear from fl4.3p that i?^(/, (?) 
is a smooth function in V^, this smoothness becomes obvious from (14.41) .) 
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Lemma 1. Let V be an open set in Mf^, and f,g ^ ^ior(^)- ^ ^ compact subset of 
V, and let eo = dist(i^, R'^ \ V) /2 and Kq =^ {a; e M"' : dist(x, K)<eQ]. Then there exists a 



constant C such that, in the notation the following estimates hold for alls G (0,eo)-' 

(i) 

(ii) 

||V/^||cO(X) < C'£°~"^||/||co."(Xo)' 

(iii) 

\\R^ {f ^ 9)\\c^\K) < C'^^°||/||co."(Rd)ll5'llco."(ifo)- 

Proof of LemmaUl Note first that, for any e G (0, ^o); -ft' is a subset of V^, so tliat g'^ and 
R'^{f,g) are well defined and smooth on an open set containing K. 

(i) For every x G we have 

- f{x)\ = I [ Q{z){f{x - ez) - fix)) dz\ 

JBi{0) 



< ll/l|co.<^(i^o) / Q{z)e''\z\'' dz 

JBi(0) 



(ii) For every x E K, we have, using the fact that g is compactly supported in -Bi(O), that 
|Vr(x)| = |-^ / Vg{^)f{y)dy\ 



^ Jv ^ 

1 

e 
1 

e 



W Q{z)f{x — ez) dz\ 

Bi{0) 

Vg{z){f{x-ez)-f{x)) dz\ 

Bi(0) 



<-\\f\\co-^iKo) [ \V0{z)\e"\zrdz 
< C£:""^||/||co.-(i^o) 

(iii) Using (14. 2 p and (14. 3 p and reasoning analogously as in the proof of (i) we obtain the 
required estimate. □ 

Lemma 2. Let V be an open set in M*^ and G Ll^^iV). We use the notation dk, 

k G {l,...,d}, to denote the partial derivative with respect to the kth variable, either in the 
classical sense, or in the sense of distributions. Suppose that 

d 

(4.5) dkfk = in the sense of distributions in V. 

k=l 
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Then, for any e > such that is non-empty, where '= {x ^ V : dist(x,]R'^ \ > e}, 
we have: 

d 

(4.6) dkfk = in the classical sense in . 

k=l 



Proof. Fix any 5 > such that is not empty. Let (f G C^(V^) be arbitrary, and observe 
that if' G C^{V), where =^ (yj * in Mf^. Assumption (14. 5 p imphes that 

d „ 

/ fk{x)dk^p%x) dx = 0. 

fc=i 

Using the fact that the regularization operator commutes with differentiation on smooth 
functions (see for instance |8l Chapter 3]), then Fubini's Theorem, and then integration by 
parts, we obtain 

d „ 

= 5^ / fk{x){dkipy{x)dx 

k=l Jy 

= ^ y fk{x) (^j dk(p{y)g%x - y) dy^ dx 

= Yj ^k^^y) {^j fkix)0'iy - ^) dx^ dy 

d p 

= Y fk{y)dk^{y) dy 



k=l 



Y^^fk (y) v{y)dy. 



.fc=i 



Since if G C^{V^) was arbitrary, the required conclusion (14. 6 p follows. 



□ 



After these preliminaries on regularization, we are now in a position to give the proof of 
our main result. 

Proof of TheoremUl We prove first the equivalence of (ii) and (iii) which, as we shall see, is 
valid for any a G (0, 1]. 

Suppose that (ii) holds. Let tp G Cl'^^(D^) be such that ([32D and ([23]) hold, where F G 
C^'"{[po,0]). Defining h as in Section 2, we then have that h G C^'^^R), and the formulae 
(I2.13p - (l2.14p are still valid. (We regard (I2.14p as a relation between the classical derivatives 
of any function with respect to the (x, y) variables and those with respect to the (g, p) 
variables, and do not assign to it any meaning in the sense of distributions.) Clearly (13.2bp - 
(13.2dp imply ( ]3.3bp -( 13.3cp . and (12. 4p implies (I2.12p . Let us now write explicitly the weak form 
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of fl3.3ap . which we need to prove: for any G C\{K)^ 

(4.7) + ^^^v ^ i^' '^'^'^^ ^ °' 

For any such (p, let (/? G Cl^D^^) be given by ip{x,y) = (p{x, —ip{x,y)) for all {x,y) G Drj. By 
changing variables in the integral, using fl2.13p - fl2.14p . one can rewrite fl4.7p as 

(4.8) / T{-ij)ify - iij^ijy)^. + i^Dfy dxdy = 0. 

But fl4.8p is valid, as a consequence of fl3.2ap . This shows that fl3.3ap holds. We have thus 
proved that (iii) holds. 

Suppose now that (iii) holds. Let h G Clf^{R) be such that ([33D and fl2:T2|) hold, where 
r G C°'°([j>oj 0]). Defining j] and t\} as in Section 2, we then have that r/ G Cp^"(]R) and 

G Ci^^(A^), and the formulae fl27[3|l - fl27[4l ) are still valid. Clearly fl3:3b|l - flOc|) imply 
f ]3.2bp -f r3.2dp . and fl2.12p implies f l2.4p . The weak form of f l3.2ap . which we need to prove, is 
written explicitly as fl4.8p . for any ip G C\{Dri)- For any such let (p G C\{K) be given by 
'Pil^p) = ¥^{(1^ h{q,p)) for all {q,p) G R. By changing variables in the integral, using fl2.13p - 
fl2.14p . one can rewrite (14. 8 p as (14.71) . But (14. 7p is valid, as a consequence of fl3.3al) . This 
shows that fl3.2ap holds. We have thus proved that (ii) holds. 

We now prove the equivalence of (i) and (ii), making essential use of the assumption a > 
1/3. 

Suppose that (i) holds. Since r] G C^'^iM) and u,v E Cp'^{Dn), it follows from (13. lap , by 
arguments similar to those in [H Lemma 3] , in which our Lemma 2 plays a key role, that there 
exists ip G Cp'^{Dr^), uniquely determined up to an additive constant, such that (12. 3p holds. 
Clearly, (12. 2 p implies (12. 4p . Also, it follows from f l3.1dp and f l3.1ep that ip is constant on each 
oi y = and y = t]{x). The additive constant in the definition of ip may be chosen so that 
f l3.2cp holds, and then f l3.2bp also holds for some constant po < 0. Using the definition of ijj 
we rewrite fl3.lbp - fl3.lcp in the weak distributional form (with ipx,'^y in the classical sense): 

(4.9a) {^l)x - {^x^y)y = -P'x in D^j, 

(4.9b) -{^Px^y)x + {^Pl)y = -Py-g in Dr,. 

Let us denote 

(4.10) F''^' P + ^\Vi>\' + gy in D,. 

It follows from (14. 9p that we have, in the sense of distributions (with ipx,'^y in the classical 
sense): 

(4.11a) Fx = ^{^l-^l)x + {i^xi^y)y mD^, 

(4.11b) Fy = {i;x^y)x-lii^l-^'y)y iuD,. 

We now show that there exists a function F G C°'°([po, 0]) such that 
(4.12) F(x,2/) =F(-^(x,y)) for all (x, 2/) G/^,. 
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Let us consider again the transformations fl2.9p -f ETTT]) . and note that fl2.13p - fl2.14p are vahd 
under the present regularity assumptions. Let F : i? — )■ M be given by F{q,p) = F{q, h{q,p)) 
in R, which is equivalent to F{x,y) = F(a;, — y)) in D^. Then our desired conclusion 
KWf is that 

(4.13) F{q,p) = T{p) for all G i?, 

for some F G C^'^'dpo, 0]). To this aim, we shall prove that, for any ip G Cq{R), 

(4.14) / Fipgdqdp = 0, 

Jr 

which, together with the condition F G C^'^{R), will imply f l4.13p for some function F G 
C°'°([po7 0]). For any such (f, let ip G Cq(D^) be given by ip{x,y) = (p{x, —ip{x,y)) for all 
(x, y) G Dr^. By changing variables in the integral, using f l2.13p and fl2.14p . fl4.14p can be 
rewritten as 



(4.15) / F {ipyifx - iix^y) dxdy = 

Thus our aim is to prove fl4.15p for any ip G Cq(Z^^). Note for later reference that this 
statement can be written in the sense of distributions (with ipx^i^y in the classical sense) as 

(4.16) {F^y)x - {F^,)y = inD^. 

Note also that, for any function 6 of class C^, a direct calculation shows the identities 



(4.17a) ^^el-9l), + {9Jy)y = 9.,A9, 

1 



(4.17b) ^ejy)x-^{9l-9l)y = 9yA9. 

(The expressions in the left-hand side of fl4.17p are similar to those occurring in the right-hand 
side of (14. lip , however (14.1 7p cannot be applied with 9 := tp, since ip is not of class in 

D^.) Let V = Drj and let, for any e > 0, let V = {{x,y) G V : dist((x, y), \V) > e}. 
Using Lemma [2] the system (14.111) implies that, for any s > such that is not empty, in 
the notation (14. ip . 

= - 0^1)% + + in^ii^.^i^.). - lR'{ijy,ijy). + R'{^.,^y)y in V\ 

F'y = m^l). - - {ryf)y + R'{^..^y). ' ^i?^(^.,^.), + {i^y.i^y)y m . 

Using the identities f l4.17p with 9 := ip^ , the above can be rewritten as 

(4.18a) F| = ^IA,P' + ^i?^(^.,^.). - ]^R'{iljy,^y)x + R'{ip:,.^y)y in V , 

(4.18b) = ryM' + i?'(^x,^,)x - ^i?^(V'x, V'x), + \R%i^y.^y)y m V . 

Let if G Cq(D^), arbitrary, and let K =^ suppy^. Let Eq ^= dist(i^, \ V)/2 and Kq =^ 
{(x, y) G : dist((a;, y),K) < Eo}. Note that is a subset of V"^, for any e G (0, Eq). Aiming 
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to prove fl4.15p for we write, for any e G {Q^Sq), 
F i^pyfx - i^x^y) dxdy 

= [ {Fijy-F'ryW^-{F^px-F'rxWydxdy+ [ F'^y^, - F'^j'^^y dxdy 
Jk Jk 

= Ie + Je. 

It is a consequence of Lemma l(i) that — )■ as £ — )■ 0. To estimate jTe, we first integrate 
by parts, then use fl4.18p to cancel some terms, and then integrate by parts again, to get 



Je 



{F^^P:, - F^ijl)^ dxdy 



K 



K 



]^R'{iIj^, '^y)x + ^^(^x, '^y)y 



+ 



K 



R'{lp^,i)y)^ - ]^R'{lp^,i)^)y + ^R%ijy,ijyjy 
1 ^ , . 1 



K 



1 

2' 



{%¥^) dxdy 
{'iplf) dxdy 



^R'{i)^,i)^) - -R'{il)y,i)y 



Wy^)x + {ijl^)y\ dxdy 



R'iipx,ipy)[i'^l^)y - i'^l^)x] dxdy. 



K 



Expanding the square brackets, we write as a sum of six terms, all of which can be estimated 
in a similar way, by using Lemma 1, to the one shown below: 



K 



R%^x,i^y){i^l^)xdx\ = \ / R'{ij^,i;y){^plip^ + ijlyLp)dx\ 



K 



where C is a constant which depends on (yj, but is independent of e G (0, Eq). The assumption 
a > 1/3 now implies that — )■ as e — )■ 0. We have thus proved that fl4.15p holds for any 
ip G Cl{Dr). As discussed earlier, this implies the existence of F G C°'°([po5 0]) such that 
fl4.12p holds. It therefore follows from (14. lip that, in the sense of distributions (with ipx,'^y 
in the classical sense), 



.3q-1 I 



(4.19a) 
(4.19b) 



^{-i^)x = ^(V'x - i^l)x + {i^xi^y)y m D 



^i-^)y = {'^x^y)x - -{i^l - ^l)y in Dr, 



Rearranging (14.19bp gives exactly (I3.2ap . Also, recalling f l4.10p . we obtain from (IS.lfp the 
validity of ( ]3.2dp for some constant Q. We have thus proved that (ii) holds. 

Suppose now that (ii) holds. We define in the velocity (m, v) by (12. 3p and, up to an 
additive constant, the pressure P, by 



(4.20) 



P = -l\Vij\'-gy + T{-^) inD,. 



Then u,v,P G C^i^{D^). Moreover, the definition of u and v implies (13. lap , while (13.2bp - 
( ]3.2dp imply ( ]3.1dp -( l3TTfj) . provided the additive constant in the definition of P is chosen 
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in a suitable way. Also, f l2.4p implies fl2.2p . It therefore remains to prove the validity of 
( I3.1b|) -( l3.1cp . Using our definition of u,v and P, ( I3.1b|) -( l3.1cp can be equivalently rewritten 
as f l4.19p . However, ( ]4.19bp is exactly fl3.2ap . which we are assuming to hold, and therefore 
it only remains to prove fl4.19ap . We now show that (14.19bp implies f l4.19ap . For notational 

convenience, we denote F =^ r(— ■?/;). We claim that, with this definition of F, fl4.16p neces- 
sarily holds. Indeed, fl4.16p can be written explicitly as f l4.15p for any ip G Co(-D^), which, 
using the same notation as earlier in the proof, is equivalent to fl4.14p for any (p G Co(-R), 

which is clearly true with our definition of F. Let V *== and let, for any e > 0, let 

=^ {{x,y) G V : dist((a;, y), \ > e}. Using Lemma |21 fl4.16p implies that, for any 
e > such that is not empty, in the notation (14. ip . 

(4.21) (F^^p,. - {F'ijDy + R'{F, ^y)^ - R%F, ^,)y = in V'. 

Let K be any compact subset of V. It follows from (12. 4p that, for all e sufficiently small, 

ipl<0 in K. 

For such values of e, (I4.2ip can be rewritten as 



(4.22) 



= ^Te^y - Te'^R'iF, i^y)^ - R^{F, i^,)y] in K. 

Wy % 



Using again Lemma [2] and (14.19bl) . we obtain, by using fl4.17bp with 9 := ■ip'^ ^ that, for any 
e > such that is not empty. 



1 



1 



(4.23) = + i?^(V'x,^,)x - -R'{i^.,i^a)y + ^R'i^y,^: 

We deduce from fl4.23p using (14.221) that, for all e sufficiently small, 

F" 



in V^. 



R%1px,i^y)x - ^R%1px,1px)y + ^R%1py,1py)y 



(4.24) 



R%F,^,)y] inK. 



We now write explicitly the weak form of fl4.19ap . which we want to prove: for any if G Cq{D^), 

1 



(4.25) 



F^x -^{ipl- i^l)^x - iipx'ipy)^y dxdy 



0. 



Let ip G Cq{Djj), arbitrary, and let K =^ suppy?. Let Eq =^ dist(ii', \ V)/2 and Kq =^ 
{{x, I/) G : dist((x, y), K) < Eq). Note that K is a subset of V^, for any e G (0, ^o)- Aiming 
to prove (I4.25P for (/?, we write, for any e G (0,eo)) 

1 



F'^x - -i^x - ^y)fx - i^xijy)^y dxdy 



{€-^l)--Mif-{ryf) 



f[F- F^]ip, - 
Jk 

F'^X - - {ryf)Vx - {rxry)Pydxdy 



^X 



xi^y) - {'iPl%)]^y dxdy 



def 
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It is a consequence of Lemma l(i) that /C^ — )■ as e — ?■ 0. Now note that fl2.4p imphes that 
there exists e e (0, eo) and 6 > such that, for all e G (0, e), 

(4.26) ipl>d inK. 

To estimate Ce, we first integrate by parts using fl4.17al) with 9 := ip'^, then use fl4.24p to 
cancel some terms, and then integrate by parts again, to obtain, for any e G (0,e), 



K 



K 



if dxdy 



+ / ^ [W{F,,Py).,- W{F,^lj,)y\^dxdy 

JK Wy 



R%1px,i^y 



K 



R%F,i;y 



R%F,i;,) 



-R%^y,^yj 



Lf ) dxdy 



—ip dxdy. 



Thus we have written £e as a sum of five terms, all of which can be estimated in a similar 
way, by using Lemma 1, to the one shown below: 



R'{ilJx.i)y)[^^] dxdy\ 



K 



R'{lpx,1py 



K 



e 

X fy 



^x + 



i^xxi^y - ^x^: 



xy 



if \ dxdy\ 



< C{e'^°'\\ipx\\co,<^(^KQ)\\''Py\\cO''^{Ko) + 



.3a-l 



11^. 



\\^y\ 



where C is a constant which depends on if, but is independent of e G (0,£), and we have 
also used (14.261) . The assumption a > 1/3 now implies that £e — ?■ as e — )■ 0. We have 
thus proved that (14.251) holds for any ip G Co(-D^), and therefore that (I4.19ap holds. This 
completes the proof that (i) holds. □ 
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